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Transport of quantum information in linear spin chains has been the subject of much theoretical 
work. Experimental studies by nuclear spin systems in sohd-state by NMR (a natural implemen- 
tation of such models) is complicated since the dipolar Hamiltonian is not solely comprised of 
nearest-neighbor XY-Heisenberg couplings. We present here a similarity transformation between 
the XY-Heisenberg Hamiltonian and the grade raising Hamiltonian, an interaction which is achiev- 
able with the collective control provided by radio-frequency pulses in NMR. Not only does this 
second Hamiltonian allows us to simulate the information transport in a spin chain, but it also 
provides a means to observe its signature experimentally. 



In many solid-state proposals for quantum computers, 
one essential task is the transport of information over 
relatively short distances. Relying on photons requires a 
frequent exchange of information between solid-state and 
light qubits that could be too costly. Linear chains of 
spins coupled by the XY-Heisenberg Hamiltonian could 
be a viable alternative P, [3] . Various schemes have been 
proposed to obtain perfect state transfer, involving en- 
gineering of the system Q or of the Hamiltonian 0, Q . 
Usually these schemes are based on nearest-neighbor cou- 
plings only, which enable analytical solutions of the dy- 
namics . Taking into account the full couplings 0, [§| 
complicates the dynamics and steps toward control se- 
quences to reduce the couplings to nearest-neighbors have 
been taken [^, [loj . A quantum simulation of this trans- 
port mechanism could connect these theoretical studies 
to the phenomenon of spin diffusion observed in dipo- 
larly coupled spin systems. While optical lattices have 
been proposed as versatile quantum simulators for a va- 
riety of spin Hamiltonians [isj , here we propose to 
use a system closer to the simulated one, taking advan- 
tage of the control techniques developed by the NMR 
community. 

In this letter we show how well-known NMR pulse se- 
quences can be used to experimentally study the trans- 
port of quantum information, in collectively controlled, 
room temperature linear spin chains. Quasi-ID spin 
systems are available in some materials like apatites 
(either hydroxyapatites or fluorine containing apatites 
l4 . llSl . Ilq). We first introduce a Hamiltonian (the so- 
called double-quantum (DQ) Hamiltonian) that is con- 
nected to the XY-Hamiltonian by a similarity transfor- 
mation and that is obtained experimentally by sequences 
of pulses and delays in dipolar spin systems. We show 
how this Hamiltonian can simulate the transport dynam- 
ics for relevant initial states, by comparing the probabil- 
ity of transfer for the two Hamiltonians. We show fur- 
thermore that the DQ-Hamiltonian enables the detection 
of successful transport, even when it is not possible to 
measure single nuclear spins. We conclude by discussing 
some of the questions that the experimental studies en- 



abled by this method could attack. 

Transport of Zeeman and dipolar energy in spin sys- 
tems, caused by energy conserving smn flips, has been a 
long term interest in spin physics [111 . 1 1711 and more re- 



cently in quantum information science |18l . Il9l . 1201 . |21 1 
The evolution of the system induced by the secular dipo- 
lar Hamiltonian 



(1) 



is coherent (as prov ed by the observation of polarization 
echoes 2^ |23|, |2J, |25|, l26|), but the complexity of the 



interaction encodes information about the many-body 
states created in obscrvablcs that are not directly measur- 
able in NMR experiments. The free-evolution transport 
presents the signature of an apparent diffusive behavior, 
with an effective decay time Tj* much shorter than spin- 
lattice relaxation times. Refocusing experiments can ex- 
tend the transverse relaxation T2 to much larger values, 
although still quite far from the theoretical limit of the 
longitudinal relaxation Ti. If no special assumptions are 
made on the values of the dipolar couplings , spin diffu- 
sion (that is, evolution under the secular dipolar Hamil- 
tonian) cannot transfer the polarization from the first to 



the last spin with any appreciable efficiency [271 . 
for a small number of spins in linear chains. On the other 
hand, a simpler Hamiltonian, the XY interaction. 



(2) 

has been studied extensively for perfect transport pur- 
poses, but it is not found naturally in spin systems. It 
has been shown that perfect transport of a state is pos- 
sible only for chains up to 3 spins 0], but if one can 
engineer the coupling strength [29| or add external ma- 
nipulation of all the spins [30[ or even just of the spins 
at the chain end [3ll |. perfect transport is achievable for 
chains of arbitrary length. 

Unfortunately, given the dipolar Hamiltonian ([T]) it is 



2 



impossible with only collective control over the spins to 
break its symmetry and obtain the XY-Hamiltonian (by 
eliminating the term cx cr^cr^). This collective control, by 
means of a sequence of radio-frequency (rf) pulses and 
delays [s^], allows us to create the DQ-Hamiltonian: 



their Fourier transform [SS, 35, 3; 



(3) 

A simple unitary transformation links the DQ- and 
XY-Hamiltonian if we restrict the couplings to nearest- 
neighbor spins only (33| . The transformation however 
cannot be obtained by an rf pulse sequence, since it is 
not a collective operation: 



dq 



(4) 



where the sum is restricted to either even or odd spins in 
the chain. 

Since the two Hamiltonians are related by a similar- 
ity transformation, one can use the DQ-Hamiltonian (ex- 
perimentally available) to simulate the dynamics of the 
transport Hamiltonian, provided one can prepare initial 
states and detect observables which arc invariant under 
the transformation. With these assumptions, even if the 
similarity transformation is not explicitly applied, 
the evolution of these particular initial states is the same. 
For thermally polarized spin systems at T > lOOmK, the 
most common initial state is to have just one spin polar- 
ized, p(0) = CT°. Since this state is invariant under the 
transformation the DQ-Hamiltonian opens the pos- 
sibility to experimentally study spin transport in dipo- 
larly coupled systems. Incidentally, the pure state usu- 
ally considered in theoretical studies, that is, the one-spin 
excitation state, is not invariant under the transforma- 
tion but since here we are considering an experimental 
implementation, mixed states are more relevant. 

Notice that the isomorphism between the XY- and DQ- 
Hamiltonian is valid also in two and three dimensions 
(and in general on hypercubes of any dimensions, which 
have already been shown to allow perfect state transport 
under the XY interaction [1^). These spin systems can 
be thought as tree- like graphs [s^l, where the first spin 
(a spin at one edge of the hypcrcube) is at the top of the 
tree, and has D neighbors in the following level, where 
D is the dimension of the hypcrcube; in general each 
node has D links to the upper level and D links to the 
lower level. The similarity transformation is performed 
by flipping all the spins of either the odd or even levels. 

We study the polarization transfer given by the XY- 
and DQ-Hamiltonian, when the transformation is not 
explicitly applied (as it will be experimentally). The dy- 
namics is solved in terms of fermion operators and of 



n K>7; flfc 
fe=i 



N 



N- 



-^sin(fcj)cj (5) 



where k = j^q^, n G integers. The nearest-neighbor 
coupling XY-Hamiltonian is diagonal in terms of fermion 
operators: 

Hxy = d^{c]cj+i+Cjc]^-^) = 2d^cos(fc)4afc (6) 



In order to diagonalize the DQ-Hamiltonian, we need a 
further transformation to Bogoliubov operators (STj : 



(7) 



where = sgn(fc). We finally obtain a diagonal form, 
T^dq = -2dJ2 cos k{dldk + dl^d^k ~ 1), (8) 

k 

with the same eigenvalues as the XY-Hamiltonian. 
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FIG. 1: Polarization transfer from spin one to spin 
N for a chain of 21 spins with nearest-neighbor coupling 
only; the results of evolution under double-quantum and 
XY-Hamiltonian are superimposed. The dipolar coupling 
strength considered are those found in a Fluorapatite crystal 
[3^ . In the inset, different behavior for even or odd number 
chain lengths 



Assuming that the polarization resides initially just on 
spin a, the initial state is: 



^ 2 

Pa{Q) = ^ - j;^— sin (ka) sin (ha) alah (9) 

k,h 

which evolves under the XY-Hamiltonian ([6]) as: 
P7it) = \- 7m '^k,h sin (^a) sin {ha) 

w „ — 2zfi(cos /c— cos h)t ^ ^ ' 
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At a time t, the polarization transferred to another spin 
b is given by Tr [p^J{t)a''^] , that is: 



(iV + 1) 



—2 I sin (fca) sin {kb)e 



(11) 



where ipk{t) = 2dtcosk. Evolution under the DQ- 
Hamiltonian with the same initial state yields: 

Pa^(*) = T - IvTT T.kM Sin (ka) sin {ha) 
isiml;k{t) cos iljh(t){aia'' - a^kah) 



«sin V'fc(i) cos ■4ih{t){a\a^_k 

' s\nipk{t) sh\^})h{t) aha\] 



-{cos iik{t) cos iph{t) a\ah 



and the 



(12) 

polarization transferred to the h spin at time t is 



(n + iy 



-Re 



sin (fca) sin {kb)e 



iil>k(t) 



(13) 

where Re[.] is the real part of a complex number. No- 
tice that P^l{t) = |l^'ab(^)ll- The polarization transfer is 
the same if the difference 6 — a is even, while there is a 
difference for — a odd [see Fig. ([T])]. In the latter case 
the observable is no longer invariant under the similarity 
transformation; to recover the result P^l{t) = -P^h (i) one 
should have applied the transformation U^^ explicitly. 

More generally, if we measure a different observable 
than cr^, wc expect different behaviors whether the evo- 
lution is driven by Ti^q or TYxy In particular the XY- 
Hamiltonian conserves the total magnetic number (the 
system remains in the zero-quantum coherence man- 
ifold), while the DQ-Hamiltonian can create multiple 
quantum coherence (MQC, states that show a coher- 
ence between different z- magnetic moment states). In 
the nearest-neighbor coupling limit, the DQ-Hamiltonian 
creates only zero and double-quantum coherences (33| 
and its dynamics can be solved exactly. (For couplings 
among all the spins Tidq creates all the even coherences, 
as observed in MQC experiments in NMR [3^V 

This property of the DQ-Hamiltonian provides a means 
to detect the occurred transfer of polarization. The 
transport of polarization cannot be detected directly (un- 
less one could introduce very strong magnetic field gra- 
dients or if one could perform single-spin detection). In 
the case of double-quantum dynamics, however, there is 
a correlation between the coherence intensities and the 
transport of the spin state, which arise from boundary 
effects in finite length spin chains 

MQC intensities can be detected in NMR, by encoding 
the coherence order into a phase associated with rotations 
around the z axis [13, S^. In MQC experiments, states 
with higher coherence orders are first prepared, for ex- 
ample using Hdq, p{0) -5" p = J2q P^''\ where q indicates 
the coherence order; after possible manipulations that de- 
pend on the goal of the actual experiment, the system has 
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FIG. 2: MQC intensities and polarization transfer in a 21- 
spin chain. The initial state was the one we can prepare 
experimentally, po = ul + [36| . 



to be brought back to a single-coherence state in order 
to be measured. The experiment is repeated Mq times, 
each time imposing a ipk = 2-Kk/Mq rotation around the 
z axis, so that different coherences acquire a phase pro- 
portional to their coherence order = qip. The signal is 
then Fourier-transformed with respect to this phase, to 
obtain the intensity of each coherence. Restricting the 
interaction to nearest-neighbor only in ID, the zero- and 
double-quantum intensities Ja{t) ~ Tr [(p*-"-*)^] can be 
calculated analytically [3^, yielding for an initial state 
with only spin a polarized : 



IF Sfe,/i sin {kaf sin {ah) 



{N + 



(14a) 



X cos [lpk{t) +-4)h{t)Y 

•^2 (0 = JnTW ^''^h sin {ka)^ sin {ha)^ 
xsiniMt) +Mt)f 

These intensities present a beating every time the po- 
larization reaches spin N + I — a and back to a. These 
beatings are particularly clear for the transfer from spin 
1 to spin N, although they would exist for the magnetiza- 
tion starting at any spin in the chain. If one is therefore 
interested in the transfer of polarization from one end of 
the chain to the other, it is possible to follow this transfer 
driven by the DQ-Hamiltonian by measuring the MQC 
intensities. If the measurement only allows to detect the 
total magnetization, the MQC intensities are then calcu- 



lated from Tr 



p{t)^">{U6qE^zU^ 



t 



resulting in 



Jo it) 



{N + 1) 



— ^sin^ (fca)cos^2V'fc(i) (15a) 



^^2"^*^ = JnTT) ^ ^^""^ ''''' ^"^'^^^ ^^^^^ 

and a less visible signature (as shown in fig. [21). 

It has been shown [3^ that with collective coherent and 
incoherent control it is possible to create the initial state 
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cr^ + (even if it is not possible to break the symmetry general problems in condensed matter physics. 



between the two spins at the end of the chain). In Fig. 
(21) we show the transport and detection method results 
for this particular initial state, realizable experimentally. 
The beatings of the MQC intensities is now faster than 
the transfer of polarization, since the coherences start 
spreading out from the two opposite ends of the chain and 
an extremum in the MQC intensity is created when the 
two waves meet at the center of the chain as well as when 
they bounce off the boundaries. This classical interfer- 
ence occurs with a positive or negative phase, depending 
on the number (odd or even respectively) of spins in the 
chain. Every two beatings, however, the maximum of 
the zero-quantum intensity correspond to the transport 
of polarization from one end to the other one. This is 
the experimentally measurable signature that transport 
of polarization has occurred. 

An experimental study of this scheme can be imple- 



mented in single crystal of fluorapatite 
These crystals present a quasi-lD structure, where flu- 
orine spins- 1/2 are arranged in linear chains, with the 
cross-chain couplings much smaller than the in-chain cou- 
plings (40:1 in a particular crystal orientation). Good 
control, long decoherence times and the availability of 
the desired initial state make this system very promising. 
In addition, a real system will provide more insight into 
the limitations of the nearest-neighbor coupling approxi- 
mation. There already exist schemes for reducinga long- 
distance interaction to nearest-neighbors only (ol llC| , but 
they are valid only for a restricted number of spins. On 
the other hand it is interesting to study the role of the 
long-distance couplings in accelerating or impeding the 
transport, as well as the cross-chain couplings. Next- 
nearest neighbor couplings and cross-chain couplings of- 
fer additional pathways that can result in an acceleration 
of information transport, which has no classical counter- 
part, as already observed in the transport of dipolar en- 
ergy in spin diffusion [l^, [2^ . It will be interesting to 
investigate the differences between the predicted rate of 
transport and the experimental one, to observe the ef- 
fects of the real couplings on the spin dynamics, as we 
move from a solvable system to systems that display fully 
quantum-mechanical dynamics. 

In conclusion, we have shown how the DQ-Hamiltonian 
simulates the transport of polarization and enables its de- 
tection. This is made possible by the observation that not 
only is this Hamiltonian related to the XY interaction via 
a similarity transformation, but it also creates coherent 
states, whose intensities are correlated to the transport of 
polarization and can be measured experimentally. With 
this scheme it will be possible to study experimentally, 
in solid-state NMR systems, the transport of polarization 
beyond exactly solvable models and explore the appear- 
ance of quantum coherence and interference effects. We 
anticipate that the control and measurement techniques 
described in this letter will allow the simulation of more 
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